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Abstract 



We extend the model of rational bubbles of [ Blanchard, 1979 ] and [ Blanchard and Watson, 1982 1 



to arbitrary dimensions d: a number d of market time series are made linearly interdependent 
via d x d stochastic coupling coefficients. We first show that the no-arbitrage condition im- 
poses that the non-diagonal impacts of any asset i on any other asset j ^ i has to vanish on 
average, i.e., must exhibit random alternative regimes of reinforcement and contrarian feed- 
backs. In contrast, the diagonal terms must be positive and equal on average to the inverse of 
the discount factor. Applying the results of renewal theory for products of random matrices to 



stochastic recurrence equations (SRE), we extend the theorem of [Lux and Sornette, 1999 1 and 
demonstrate that the tails of the unconditional distributions associated with such rf-dimensional 
bubble processes follow power laws (i.e., exhibit hyperbolic decline), with the same asymp- 
totic tail exponent /i < 1 for all assets. The distribution of price differences and of returns is 
dominated by the same power-law over an extended range of large returns. This small value 
/i < 1 of the tail exponent has far-reaching consequences in the non-existence of the means 
and variances. Although power-law tails are a pervasive feature of empirical data, the numer- 
ical value fi < 1 is in disagreement with the usual empirical estimates /i w 3. It, therefore, 
appears that generalizing the model of rational bubbles to arbitrary dimensions does not allow 
us to reconcile the model with these stylized facts of financial data. The non-stationary growth 



rational bubble model seems at present the only viable solution [Sornette, 2001 1 



1 The model of rational bubbles 



[Blanchard, 197S] and [Blanchard and Watson, 1982] originally introduced the model of rational 
expectations (RE) bubbles to account for the possibility, often discussed in the empirical literature 
and by practitioners, that observed prices may deviate significantly and over extended time intervals 
from fundamental prices. While allowing for deviations from fundamental prices, rational bubbles 
keep a fundamental anchor point of economic modelling, namely that bubbles must obey the con- 
dition of rational expectations. In contrast, recent works stress that investors are not fully rational, 
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or have at most bound rationality, and that behavioral and psychological mechanisms, such as herd- 
ing, may be important in the shaping of market prices [Thaler, 1993, Shefrin, 2000 , [Shieifer, 200C ] . 
However, for fluid assets, dynamic investment strategies rarely perform over simple buy-and-hold 
strategies [ Malkiel, 1999| ], in other words, the market is not far from being efficient and little arbi- 
trage opportunities exist as a result of the constant search for gains by sophisticated investors. Here, 
we shall work within the conditions of rational expectations and of no-arbitrage condition, taken 
as useful approximations. Indeed, the rationality of both expectations and behavior often does not 
imply that the price of an asset be equal to its fundamental value. In other words, there can be 
rational deviations of the price from this value, called rational bubbles. A rational bubble can arise 
when the actual market price depends positively on its own expected rate of change, as sometimes 



occurs in asset markets, which is the mechanism underlying the models of [Blanchard, 1979] and 



[Blanchard and Watson, 1982] 



In order to avoid the unrealistic picture of ever-increasing deviations from fundamental values, 



[ [Blanchard and Watson, 1982| ] proposed a model with periodically collapsing bubbles in which the 
bubble component of the price follows an exponential explosive path (the price being multiplied by 
at = a > 1) with probability tt and collapses to zero (the price being multiplied by a t = 0) with 
probability 1 — n. It is clear that, in this model, a bubble has an exponential distribution of lifetimes 
with a finite average lifetime vr/(l — it). Bubbles are thus transient phenomena. The condition of 
rational expectations imposes that a = 1/5, where 5 is the inverse of the discount factor. In order 
to allow for the start of new bubbles after the collapse, a stochastic zero mean normally distributed 
component b t is added to the systematic part of B t . This leads to the following dynamical equation 



+i 



a t X t + b t , 



(1) 



where, as we said, at = a with probability it and at = with probability 1 — it. There is a huge 
literature on theoretical refinements of this model and on the empirical detectability of RE bubbles 
in financial data (see [Camerer, 1989] and [ Adam and Szafarz, 1992| ], for surveys of this literature). 

Recently, [ f^ux and Sornette, 1999 ] studied the implications of the RE bubble models for the un- 
conditional distribution of prices, price changes and returns resulting from a more general discrete- 
time formulation extending ([I]) by allowing the multiplicative factor at to take arbitrary values 
and be i.i.d. random variables drawn from some non-degenerate probability density function (pdf) 
P a {o)- The model can also be generalized by considering non-normal realizations of b t with dis- 
tribution Pb{b) with E[b t ] = 0, where E[-] is the expectation operator. Since in ([!]) the bubble X t 
denotes the difference between the observed price and the fundamental price, the "bubble" regimes 
refer to the cases when X t explodes exponentially under the action of successive multiplications 
by factor at, at+i, ■■■ with a majority of them larger than 1 but different, thus adding a stochastic 



component to the standard model of [ [Blanchard and Watson, 1982| ]. 

Provided E[ln a] < (stationarity condition) and if there is a number \i such that < Eflb^] < 
+oo, such that 



E[|aH = 1 



(2) 



and such that E[[a| M In |a|] < +oo, then the tail of the distribution of B is asymptotically (for large 
X's) a power law QKesten, 1973] poldie, 199 1| ] 

Px{X) dX ( 



\X\ l +v 



dX 



(3) 



with an exponent \i given by the real positive solution of (^). Rational expectations require in 
addition that the bubble component of asset prices obeys the "no free-lunch" condition 



5 • E[X 



t+U 



X t 



(4) 



where 5 is the discount factor < 1. Condition (Q) with (|l|) imposes the condition 

E[a] = 1/6 > 1 , (5) 

on the distribution of the multiplicative factors a t . Since the function E[|a| M ] is convex, [ Lux and Sornette, 1999 ] 
showed that this automatically enforces \i < 1. It is easy to show that the distribution of price differ- 
ences has the same power law tail with the exponent n < 1 and the distribution of returns is domi- 
nated by the same power-law over an extended range of large returns [ Lux and Sornette, 1999 ]. Al- 
though power-law tails are a pervasive feature of empirical data, these characterizations are in strong 
disagreement with the usual empirical estimates which find \i m 3 [ de Vries, 1994 , Lux, 1996 



Pagan, 3uillaume et al., 1997, Gopikrishnan et al., 1998]. [ Lux and Sornette, 1999j ] concluded that 
exogenous rational bubbles are thus hardly reconcilable with some of the stylized facts of financial 
data at a very elementary level. 



2 Generalization of rational bubbles to arbitrary dimensions 
2.1 Introduction 

In reality, there is no such thing as an isolated asset. Stock markets exhibit a variety of inter- 
dependences, based in part on the mutual influences between the USA, European and Japanese 
markets. In addition, individual stocks may be sensitive to the behavior of the specific industry 
as a whole to which they belong and to a few other indicators, such as the main indices, interest 
rates and so on. [ Mantegna, 1999j , Bonanno, Lillo and Mantegna, 2001 ] have indeed shown the ex- 



istence of a hierarchical organization of stock interdepences. Furthermore, bubbles often appear 



to be not isolated features of a set of markets. For instance, [ [Flood et al., 1984| ] tested whether a 
bubble simultaneously existed across the nations, such as Germany, Poland, and Hungary, that ex- 
perienced hyperinflation in the early 1920s. Coordinated bubbles can sometimes be detected. One 
of the most prominent example is found in the market appreciations observes in many of the world 



markets prior to the world market crash in Oct. 1987 [Barro et al., 1989]. Similar intermittent coor- 



dination of bubbles have been detected among the significant bubbles followed by large crashes or 



severe corrections in Latin- American and Asian stock markets [lohansen and Sornette, 2000]. It is 
therefore desirable to generalize the one-dimensional RE bubble model ([[]) to the multi-dimensional 
case. One could also hope a priori that this generalization would modify the result /i < 1 obtained 
in the one-dimensional case and allow for a better adequation with empirical results. As we shall 
show, this turns out to be ill-born for reasons that we shall explain in details. 

The simplest such generalization is to consider two assets X and Y with prices respectively 
equal to X t and Y t at time t, evolving according to 

X t+1 = a t X t + b t Y t + m (6) 

Y t+1 = c t X t + d t Y k + e t (7) 

where a t , b t , ct and dt are drawn from some multivariate probability density function. The two 
additive noises rjt and et are also drawn from some distribution function with zero mean. The 
diagonal case bt = ct = for all t recovers the previous one-dimensional case with two uncoupled 
bubbles, provided rj t and et are independent. 

Rational expectations require that X t and Y t obey both the "no-free lunch" condition (f|), i.e., 
S • EfXt+i] = X t and 5 • E[Y t+1 ] = Y t . With this gives 

(E[at] - r 1 ) X t + E[bt]Y t = 0, (8) 
E[c t ]X t + (E[d t ]-5- 1 )Y t = 0, (9) 



where we have used E[r/ t ] = E[e t ] = 0. The two equations (|||9|) must be true for all times, 
i.e. for all values of X t and Y t visited by the dynamics. This imposes E[b t ] = E[c t ] = and 
E[at] = E[dt] = 8~ l . We are going to retrieve this result more formally in the general case. 

2.2 General formulation 

A generalization to arbitrary dimensions lead to the following stochastic random equation (SRE) 

X t = A t X t _i + B t (10) 

where (X t , B t ) are d-dimensional vectors. Each component of X t can be thought of as the price 
of an asset above its fundamental price. The matrices (A t ) are identically independent distributed 
d x d-dimensional stochastic matrices. We assume that B t are identically independent distributed 
random vectors and that (X t ) is a causal stationnary solution of (|T0|). Generalizations introducing 
additional arbitrary linear terms at larger time lags such as Xt-2, • •• can be treated with slight mod- 
ifications of our approach and yield the same conclusions. We shall thus confine our demonstration 
on the SRE of order 1, keeping in mind that our results apply analogously to arbitrary orders of 
regressions. 

To formalize the SRE in a rigorous manner, we introduce in a standard way the probability 
space (fi, J 7 , P) and a nitration (Ft)- Here P represents the product measure P = Fx <8> Fa <8> P_b> 
where Fx, Fa and P#, are the probability measures associated with {X t }, {A t } and {B t }. We 
further assume as is customory that the stochastic process (X t ) is adapted to the filtration (Tt). 

In the following, we denote by | • | the Euclidean norm and by 1 1 • 1 1 the corresponding norm for 
any d x <i-matrix A 

||A|| = sup |Ax|. (11) 

|x|=l 

In the next section ||, we formalize the "no-free lunch" condition for the SRE ([K]) and show 
that its entails in particular that the spectral radius (larger eigenvalue) of E[A t ] must be equal to the 
inverse of the discount factor, hence it must be larger than 1. In section ||, we synthetize the main 
results on the renewal theory of SRE of the type ( JTo| ) and show that the condition imposing that the 
exponent of the power law tails be larger than 1 implies that the spectral radius of E[A t ] must be 
smaller than 1. By the reverse logic, the "no-free lunch" condition automatically enforces our main 
result that the exponent is less than 1. 

3 The no-free lunch condition 

3.1 No-free lunch condition under the risk-neutral probability measure 

The "no-free lunch" condition is equivalent to the existence of a probability measure Q equivalent 
to P such that, for all self-financing portfolio P t , ip^ is a Q-martingale, where So,t = Yl^Zo <^ _1 > 
5i = (1 + is the discount factor for period i and is the corresponding risk-free interest rate. 

It is natural to assume that, for a given period i, the discount rates are the same for all 
assets. In frictionless markets, a deviation for this hypothesis would lead to arbitrage opportunities. 
Furthermore, since the sequence of matrices {A t } is i.i.d. and therefore stationnary, this implies 
that 5t or rt must be constant and equal respectively to 5 and r. 

Under those conditions, we have the following proposition : 



A 



Proposition 1 

The stochastic process 

X t = A t X t _! + B t (12) 
satisfies the no-arbitrage condition if and only if 



E Q [A] = h d . (13) 
d 



The proof is given in the Appendix A. 



The condition (13) imposes some stringent constraints on admissible matrices A t . Indeed, 
while A t are not diagonal in general, their average must be diagonal. This implies that the off- 
diagonal terms of the matrices A t must take negative values, sufficiently often so that their averages 
vanish. The off-diagonal coefficients quantify the influence of other bubbles on a given one. The 



condition ( |13| ) thus means that the average effect of other bubbles on any given one must vanish. 
It is straightforward to check that, in this linear framework, this implies an absence of correlation 
between the different bubble components E[X( k )X^)] = for any k / I. 

In constrast, the diagonal elements of A t must be positive in majority in order for £f>[Ajj] = 
<JW ( for all i's, to hold true. In fact, on economic grounds, we can exclude the cases where the 
diagonal elements take negative values. Indeed, a negative value of An at a given time t would 
imply that xf ] at 
financial process. 



(i) 

imply that X t abruptly change sign between t — 1 and t, what does not seem to be a reasonable 



3.2 The no-free lunch condition under historic probability measure 

The historical P and risk-neutral Q probability measures are equivalent. This means that there exists 
a non-negative matrix h(#) = such that, for each element indexed by i, j, we have 

Ep[Aij] = EQ[hij ■ Aij] (14) 
= hij{e%) ■ EgiAij] for some 0g € R . (15) 

The second equation comes from the well known result : 

J f{6) ■ g{6) dfi(6) = g(e ) • J f{9) du.{9) for some 6 E R . (16) 



We thus get 

E^Aij] = if i ^ j (17) 

Ev[Aij] = ^ if i = j , (18) 

where the tfW can be different. We can thus write 

£ P [A]=<Fi, where F T = diag[5W" 1 ,---,^" 1 ] . (19) 

Appendix B gives a proof showing that #W 1 i s indeed the genuine discount factor for the i-th 
bubble component. 



4 Renewal theory for products of random matrices 

In the following, we will consider that the random d x d matrices At are invertible matrices with 
real entries. We will denote by GL^(R) the group of these matrices. 



4.1 Definitions 
Definition 1 (Feasible matrix) 

A matrix M G GL d (R) is V -feasible if there exists an n G N and Mi, • • • , M n G supp(P) such that 
M = Mi • • • M n and if M has a simple real eigenvalue q(M) which, in modulus, exceeds all other 
eigenvalue of M. 

Definition 2 

For any matrix M G GLd(R) and M' its transpose, MM' is a symmetric positive definite matrix. 
We define A(M) the square root of the smallest eigenvalue of MM'. 

4.2 Theorem 



We extend the theorem 2.7 of Davis et al. [Davis et al., 1999], which synthetized Kesten's theorems 



3 et 4 in [Kesten, 1973], to the case of real valued matrices. The proof of this theorem is given in 



[Le Page, 1983]. 



Theorem 1 

Let (A n ) be an i.i.d. sequence of matrices in GLd(M) satisfying the following set of conditions: 
HI : for some e > 0, E ¥a [||A|| e ] < 1 , 

H2 : For every open U C Sd-i (the unit sphere in M. d ) and for all x G Sd-i there exists an n such 
that 

PtI tt^— 4^ eu) >0. (20) 
U|xAi---A n || J 

H3 : The group {In \q(M) \,M is ¥ A -feasible} is dense in E . 

H4 : for all r G M d , Pr{Air + B x = r} < 1 . 

H5 : There exists a kq > such that 

E Pa ([A(A!)r) > 1 . (21) 

H6 : With the same kq > as for the previous condition, there exists a real number u > such 
that 

£ PA ([sup{||A 1 ||,||B 1 ||}n +u ) <oo, 
&,(||A 1 |h«)<oo. 



Provided that these conditions hold, 

• there exists a unique solution k\ G (0, kq] to the equation 

lim -ln^[||A 1 ---A n || K1 ] = 0, (23) 

n— too ti 

• If (X n ) is the stationary solution to the SRE in (jJ^) then X is regularly varying with index 
k\. In other words, the tail of the marginal distribution of each of the components of the 
vector X is asymptotically a power law with exponent k\. 



4.3 Comments on the theorem 



4.3.1 Intuitive meaning of the hypotheses 

A suitable property for an ecomonic model is the stationnarity, i.e. the solution X t of the SRE ( JlOj ) 
does not blow up. This condition is ensured by the hypothesis (HI). Indeed, E^ A [In 1 1 A| |] < im- 
plies that the Lyapunov exponent of the sequence {A n } of i.i.d. matrices is negative [ pavis et al., 1999 ]. 
And it is well known, that the negativity of the Lyapunov exponent is a sufficient condition for the 
existence of a stationnary solution X t . 

However, this condition can lead to a too fast decay of the tail of the distibution of {X}. This 
phenomenon is prevented by (H5) which means intuitively that the multiplicative factors given by 
the elements of A t sometimes produce an amplification of X t . In the one-dimensional bubble case, 
this condition reduces to the simple rule that at must sometimes be larger than 1 so that a bubble 
can develop. Otherwise, the power law tail will be replaced by an exponential tail. 

So, (HI) and (H5) keep the balance between two opposite objectives : to obtain a stationnary 
solution and to observe a fat tailed distribution for the process (X t ). 

Another desirable property for the model is the ergodicity : we expect the price process (X t ) to 
explore the entire space Mr. This is ensured by hypothesis (H2) and (H4) : hypothesis (H2) allows 
X t to visit the neighborhood of any point in M. d , and (H4) forbids the trajectory to be trapped at 
some points. 

Hypothesis (H3) and (H6) are more technical ones. The hypothesis (H6) simply ensures that 
the tails of the distributions of A t and B t are thinner than the tail created by the SRE (JTo|), so that 
the observed tail index is really due to the dynamics of the system and not to an ill-posed problem 
where the tail distribution of A t or B t is fat enought to dominate the large deviations behavior of 
the process. The hypothesis (H3) expresses some kind of aperiodicity condition. 



4.3.2 Intuitive meaning of ( 23 ) 



The equation (|23J) determining the tail exponent k\ reduces to (|2|) in the one-dimensional case, 
which is simple to handle. In the multi-dimensional case, the novel feature is the non-diagonal 
nature of the multiplication of matrices which does not allow in general for an explicit equation 
similar to (0). 



4.4 Constraint on the tail index 

The first conclusion of the theorem above shows that the tail index k\ of the process (X t ) is driven 
by the behavior of the matrices (A t ). We will then state a proposition in which we give a majoration 
of the tail index. 

Proposition 2 

A necessary condition to have k\ > 1 is that the spectral radius of Ep A [A] be smaller than 1 : 

k 1 >1=>p(E Pa [A])<1. (24) 



The proof is given in the Appendix C. 

This proposition, put together with Proposition 1 above, will allow us to derive our main result. 



5 Consequences for rational expectation bubbles 



We have seen in section || from Proposition 1 that, as a result of the no-arbitrage condition, the 
spectral radius of the matrix £p[A] is greater than 1. As a consequence, by application of the 
converse of Proposition 2, we find that the tail index k\ of the distribution of (X) is smaller than 1. 
This result does not rely on the diagonal property of the matrices E[A t ] but only on the value of the 
spectral radius imposed by the no-arbitrage condition. 



This result generalizes to arbitrary ci-dimensional processes the result of [ Lux and Sornette, 1 
As a consequence, d-dimensional rational expectation bubbles linking several assets suffer from the 
same discrepancy compared to empirical data as the one-dimensional bubbles. It would therefore 
appear that exogenous rational bubbles are hardly reconcilable with some of the most fundamental 
stylized facts of financial data at a very elementary level. 

A possible remedy has recently been suggested [Sornette, 2001] which involves an average 
exponential growth of the fundamental price at some return rate rj > 0. With the condition that 
the price fluctuations associated with bubbles must on average grow with the mean market return 
rj, it can be shown that the exponent of the power law tail of the returns is no more bounded by 1 
as soon as r f is larger than the discount rate r and can take essentially arbitrary values. It would be 
interesting to investigate the interplay between inter-dependence between several bubbles and this 
exponential growth model. 



APPENDIX A: proof of Proposition 1 on the no-arbitrage condition 

Let Pt be the value at time t of any self-financing portfolio : 

P t = Wt'Xt, (25) 

where W t ' = (W±, W d ) is the vector whose componants are the weight of the different assets 
and the prime denotes the transpose. The no-free lunch condition reads : 

Pt = 5 ■ E q [P t+1 \F t ] V{P t } t > . (26) 

Therefore, for all self-financing strategies (Wt), we have 

W t+ i / | J ES Q [A]-il d |x t = VX t GM d , (27) 

where we have used the fact that (W t +i) is {T t ) -measurable and that the sequence of matrices 
{A t } isi.i.d. 

The strategy Wt' = (0, ■ • • , 0, 1, 0, ■ • • , 0) (1 in i th position), for all t, is self-financing and 
implies 

{<Hi,<h2, • • • , an - -, • • • , a id ) ■ {xj^^xj®, • • • , xf\- ■ ■ , X^)' = , (28) 
for all X t € M. d . We have called the (i,j) th coefficient of the matrix Eq[A]. As a consequence, 

(an,a i2 , ■ ■ ■ ,au - ^, ■ ■ ■ ,a id ) = Vi , (29) 

and 

Eq[A] = ^Id . (30) 

The no-arbitrage contidition thus implies : £q[A] = jl^. 

We now show that the converse is true, namely that if Eq [A] = ll^ is true, then the no-arbitrage 
condition is verified. Let us thus assume that £q[A] = ^Ij hold. Then 

E®[P t+1 \f t ] = ^[Wt+i'-Xt+il^t] (31) 

= Wt+i'-^QpCt+il^t] (32) 

= Wt+i'-^QtAt+iXt + Bt+il^i] (33) 

= W t+1 ' ■ E Q [A t+1 \F t ] -X t (34) 

= W t+1 '-E Q [A] -X t (35) 

= ^Wt+i'-Xt. (36) 

The condition that the portfolio is self-financing is W{. +1 X t = W{.X t , which means that the 
weights can be rebalanced a priori arbitrarily between the assets with the constraint that the total 
wealth at the same time remains constant. We can thus write thus 

E Q [P t+1 \F t ] = Jwt+i'-Xt (37) 

= \Pt- (38) 

Therefore, the discounted process {Pt} is a Q-martingale. 



APPENDIX B: proof that tfW 1 is the discount factor for the i-th bubble 
component in the historical space 

Here, we express the no-free lunch condition in the historical space (or real space). The condition 
we will obtain is the so-called "Rational Expectation Condition", which is a little bit less general 
than the condition detailed in the previous appendix A. 

Given the prices {xjp}f.<t of an asset, labeled by i, until the date t, the best estimation of its 
price at t + 1 is E^[Xf\ |JF t ]. So, the RE condition leads to 



E v [X® x \F t ]-X, 



(0 



(i) 



t > 



(39) 



(i) 

where r\ ' is the return of the asset i between t and i + 1. As previously, we will assume in what 



follows that r 
i reads 



(0 



is time independent. Thus, the rational expectation condition for the assets 



X 



6® ■ E f 
5® ■ E T 



v(') I T 



(40) 
(41) 



where <5W is the discount factor. 

A priori, each asset has a different return. Thus, introducing the vector X t whose i th componant 



x 



(i) 



is -gfc-, we can summarize the rational expectation condition as 



X t — Ef [Xt+il-Ft] . 



(42) 



Again, we evaluate the conditional expectation of (|10|), and using the fact that {A t } are i.i.d., 
we have 

£k[Xt|.Ft-i]=£p[A]Xt-i. (43) 



This equation together with (P 



which can be rewritten 



leads to 

X t _i = £r>[A]Xt_i, 

E P [A] - 5^) X t _! = 



(44) 



(45) 



where S^ 1 = disLg[S^ 1 ■ ■ ■ 5^ 1 ] is the matrix whose i th diagonal component is #W 1 and 
elsewhere. 

The equation (45) must be true for every X t -i £ M. d , thus 



E F \A] 



(46) 



which is the result announced in section (3.2). 



1 1\ 



APPENDIX C: proof of Proposition 2 on the condition k\ < 1 

First step : Behavior of the function 



/(«)= lim i 1ii.E3p a [|| A*-- -Ail I" 1 ] (47) 



in the interval [0, kq]. 



In [ ]Kesten, 1973| ], Kesten shows that the function / has the following properties : 

• / is continuous on [0, kq], 

• /(0) = and /(/so) > 0, 

• /'(0) < (this results from the stationarity condition), 

• / is convex on [0, kq\. 

Thus, there is a unique solution k\ in (0, kq) such that = 0. In order to have n\ > 1, it is 

necessary that /(l) < 0, or using the definition of /: 

lim -ln£p,||A n --- AJI < 0. (48) 

n^oo n 

The qualitative shape of the function /(«) is shown in figure |l|. 

Second step : 

The operator 1 1 • 1 1 is convex: 

Va€[0,l] and V(A, C) d x d-matrices, ||aA + (1 - a)C|| < a||A|| + (1 - a)||C|| . (49) 

Thanks to Jensen's inequality, we have 

^HAn-.-AiH > H^^-.-A!]!! . (50) 

The matrices (A n ) being i.i.d., we obtain 

^HAn-.-AiH > ||{7i; PA [A]} n || . (51) 

Now, let consider the normalized-eigenvector x max associated with the largest eigenvalue 

A max = p(£ PA [A]) , (52) 

where p(Ep A [A]) is the spectral radius of A. By definition, 

|| {E FA [A]} n || > I {^[Alfx^l = A max . (53) 

Then 

lim -m£ PA ||A n ---Ai|| > lim - In p(E P A [A]) n = In P (E Pa [A}) . (54) 

n— >oo n n^oo n 

Now, suppose that p(Ep A [A]) > 1. We obtain 

/(1)= lim -ln^l^-.-AxH >0, (55) 

n— >oo n 

which is in contradiction with the necessary condition (|^). 
Thus, 

f{l)<0=>p(Er A [A])<l. (56) 
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Figure 1: Schematic shape of the function f{n) defined in (07|). 



